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In this paper we study tiie BCS Bose Einstein condensation (BEC) crossover scenario within the supercon- 
ducting state, using a T-matrix approach which yields the ground state proposed by Leggett. Here we extend 
this ground state analysis to finite temperatures T and interpret the resulting physics. We find two types of 
bosonic-like excitations of the system: long lived, incoherent pair excitations and collective modes of the su- 
perconducting order parameter, which have different dynamics. Using a gauge invariant formalism, this paper 
addresses their contrasting behavior as a function of T and superconducting coupling constant g. At a more 
physical level, our paper emphasizes how, at finite T, BCS-BEC approaches introduce an important parameter 
Apg = A^ — Age into the description of superconductivity. This parameter is governed by the pair excita- 
tions and is associated with particle-hole asymmetry effects which are important for sufficiently large g. In the 
P3 . fermionic regime, Apg represents the difference between the square of the excitation gap A^ and that of the 

O ' superconducting order parameter A^^. The parameter A^g, which is necessarily zero in the BCS (mean field) 

^ , limit increases monotonically with the strength of the attractive interaction g. It follows that there is a significant 

$_H ' physical distinction between this BCS-BEC crossover approach (in which g is the essential variable which de- 

C^' termines Apg) and the widely discussed (Coulomb-modulated) phase fluctuation scenario in which the plasma 

3 , frequency is the tuning parameter. Finally, we emphasize that in the strong coupling limit, there are important 

differences between the composite bosons which arise in crossover theories, and the usual bosons of the (inter- 
»3 ' acting) Bose liquid. Because of constraints imposed on the fermionic excitation gap and chemical potential, in 

crossover theories, the fermionic degrees of freedom can never be fully removed from consideration. 

rJ^ • PACS numbers: 74.20.-z, 74.20.Fg, 74.25.Nf 
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I. INTRODUCTION 

J> . The observation of an excitation gap above Tc (called the "pseudogap") in the underdoped cuprate superconductors has been 

^^ the focus of much current research. Presumably an understanding of this state will help unravel the formal machinery, if not the 
attractive pairing mechanism behind high temperature superconductivity. It is now widely believed that this state is associated 
with the underlying superconducting phase, in large part because the (d-wave) symmetry of the pseudogap is foundEl^El to be 
the same as that of the excitation gap and order, parameter in the T < Tc state. Among-viable candidates for the origin of the 
pseudosan-state are phase fluctuation scenariosH'Q, d-wave nodal excitation mechanismsO and a BCS Bose-Einstein "crossover 
Qs^ _ picture'Tfo. 
"Y^ ■ The goal of the present paper is to discuss the last of these, the crossover scenario within the superconducting state. Our work 

Cd ' is directed towards the fundamental issues of the crossover problem, with lesser emphasis on the physics of the cuprates. We 

H present a generalized overview based on a finite temperature T matpx formulation. Our aim is to provide a useful understanding 
i_^ and to extend the physics of the well characterized ground statellj. In the process we establish a clear distinction between 

incoherent, finite center of mass momentum, pair excitations and order parameter fluctuations (i.e., collective modes), and their 
respective dynamics. We formulate a gauge invariant description of the electrodynamic response with an emphasis on particle- 

Q hole asymmetry which is necessarily very important. Although we introduce a generalized T matjds jifuprnach, special attention 
will be paid to one particular version which has been extensively discussed in our previous workEloiljIi^. 

In the BCS Bose-Einstein condensation (BEC) crossover approach, it is presumed that there is a smooth evolution, with in- 
creasing attractive coupling constant g, from BCS superconductivity, in which strongly overlapping Cooper pairs form and Bose 

;_ j condense at precisely the same temperature Tc, to a quasi-ideal Bose gas state (to be characterized in more detail below) in which 
tightly bound fermion pairs (composite bosons) form at temperatures much higher than their Bose condensation temperature Tc- 
In this latter case theie is an excitation gap (pseudogap) for fermionic excitations well above Tc- These ideas date back to the late 
sixties when EaglesliZi first drew attention to the possibility of smoothly interpolating between tjhe BCS and Bose-Einstein ground 
state descriptions of superconductivity. This was followed by a well known paper by Leggettll3 who presented an interpolation 
scheme for the ground state based on a variational wave-function, in which the chepaical potential ^ was self consistently varied 
(with increasing g) from Ep to large negative values. Nozieres and Schmitt-Rinkli3, hereafter referred to as NSR, generalized 
these previous-approaches by presenting a crossover theory for computing Tc- Although this leading order theory was not fully 
self consistentli3, nevertheless, much of the essential physics pertaining to finite T was summarized in Ref. p^ 
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With the discovery of the short coherence length cuprates, several groups noted the relevance of this body of theoretical 
work. Indeed, this recognition was made well in advance of a community-wide appreciation of pseudogapjphenomena, which 
phenomena have only served to re-enforce interest in these xrossover schemes. Randeria and co-workersc^ were among the 
first to apply the NSR approach to the cuprates. Micpas et aK3 presented detailed studies of the attractive Hubbard model with 
varying on-site Coulomb interaction [/, and Uemuraa noted, on the basis of unusual correlations deduced from /iSR experiments, 
that the cuprates exhibited aspects of bosonic character, as might be expected in a crossover theory. Our own groupO has also 
addressed cuprate issues in the past year using the formalism of the present paper. 

Attempts to go beyond the NSR schaop-at finite T are relatively more recent and almost exclusiiiehLrestricted to twa dimen- 
sional systems. Numerical simulationsEja on the attractive Hubbard model, along with numericalEjEa and analyticalEj studies 
of the so-called FLEX schemeo have provided some insights—'Diis-diagrammatic FLEX approach should be contrasted with 
an alternative (called the "pairing appmximation") which weOEptHES have introduced into the literature and which is based on 
earlier work by Kadanoff and Martinet, and extended by PattoncZl. In contrast to the FLEX scheme this latter approach precisely 
yields BCS theo«j4ri the small g limit. What is more important in distinguishing our work from that of others, however, is 
our direct focust2Hlil on crossover effects within the superconducting state. This is the topic of the present paper, as well, and 
necessarily requires studies of systems in higher than two dimensions, where Tc is non-zero. The ultimately decisive factor in 
determining whether these crossover theories or any other alternatives are apprporiate for the cuprates, may well come from the 
predicted behavior below Tc. This fact has also been emphasized by DeutscheiEa. 

II. OVERVIEW OF BCS-BEC CROSSOVER THEORIES 

A. Physical Picture at T = 0: Previous Work 

All crossover theories are based on an underlying "generic" Hamiltonian describing the attractive interaction between fermions 
with opposite spin orientation (for spin singlet pairing) 

k<T 

+ 2^ ^k,k'Ck+q/2tC-k+q/2i^-k'+q/2J.Ck'+q/2T' (1) 

kk'q 

where cj^^ creates a particle in the momentum state k with spin a, and ek is the energy dispersion measured from the chemical 
potential /i (we take h = ks = 1). For the jellium case ek = k'^/2TO — /i, while for the lattice case we consider an anisotropic 
tight-binding model with ek = 2t|| (2— cos fca;— cos fcj,)+2f^ (1— cos /c^)— /i, where t|| (t±) is the hopping integral for the in-plane 
(out-of-plane) motion, and we set the lattice constant a = 1. Here we assume a separable pairing interaction Vk.k' = g'Pk'P'k', 
where g — — \g\ is the coupling strength and the momentum dependence of the function (p-k^, pilects the pairing anisotropy, and 
its form depends on the particular model under consideration. For the jellium case we takeEEl ipi^ — {1 + fc^/Zcp) ^^^, where 
l/fco gives the range of the interaction and represents a soft cutoff in momentum space for the pairing interaction. In the lattice 
case, (/'k = 1 for s-wave pairing symmetry, and ipi^ = cos k^ — cos ky for d-wave symmetry. To simplify the notation, until 
Sec. V, when the effects of (^k become relevant, we choose to write our equations with i/Jk = 1- 

It is important to stress the key assumptions of these crossover theories: (i) only two body fermionic interactions are included, 
(ii) In calculations of equilibrium properties (as distinct from studies of the collective mode of the superconducting order param- 
eter), repulsive Coulomb interactions between fermions are absorbed into the effective pairing interaction Vk,k' • These Coulomb 
effects are presumed to be weak enough so that the attractive interactions driving superconductivity dominate, and it is assumed 
that Coulomb interactions do not occur between fermion pairs or composite bosons. 

Using this Hamiltonian, LeggettEj found that the BCS ground state wavefunction 

I*) = Hk (uk + t'k4|cLkx) |vac) (2) 



is appropriate to both the weakly and strongly interacting limits, provided ^k = ^(1 ~ ek/^-k), u^ = i(l + Ck/i^k), where 

i and the implicit parameter ji follow from the following two self consistent 

-^ + E^-0, (3) 



-Ek = ■\/e| + A^. Moreover, the explicit parameter A and the implicit parameter ji follow from the following two self consistent 
conditions 
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and 



n = 2Y,vl, (4) 

k 

where n is the electron number density. 

In the weak coupling limit, the above equations yield the usual BCS state. By contrast, in the strong coupling (i.e., g — > oo) 
limit, for the case of jellium, the system corresponds to a Rose condensation of nonoverlapping and nourinteracting tightly 
boundr^irs of fermions, which resemble to diatomic molecules. This is an "essentially ideal Bose gas"t3. Indeed, it can be 
shownEJ that in the strong coupling limit, Eq. (3) can be recast in the form of the Schroedinger equation, written in momentum 
space, for an isolated diatomic molecule, .Gonsisting of the two fermions. The strong coupling limit was argued to be slightly 
more complicated for a lattice of fermionsta which, at moderate densities, is to be distinguished from jellium. Here a canonical 
transformation can be used to partially integrate out the fermions and the Pauli principle then leads-|to hard core repulsions 
between composite bosons. Nevertheless, even for a lattice, this scheme is viewed as reasonabletS so that the BCS wave 
function assumption, or equivalently Eqs. (Ij) and (Q), represent an effective mean field approximation to the solution of the hard 
core composite boson problem. This essentially ideal Bose gas treatment of the ground state in the strong coupling limit is to 
be contrasted with the behavior found in the collective mode spectrum, where calculations lead to results similar to those for a 
weakly interacting Bose gas. These calculations suggest the presence of an effective boson-boson interaction in this limit. 

Within crossover theories, the T = behavior of the Anderson Bogoliubov (AB) mode has been studied in Refs. ^, ^ 
and pi In the long wavelength (small q) limit the dispersion relation of the sound-like AB mode is linear: uj = cq, where c 
is the AB mode velocity. In weak coupling, the usual BCS limit is obtained for the AB mode velocity c = vp/^. In strong 
coupling, c reflects an effective boson-boson interactional which arises from the Pauli principle. We caution that even though the 
phrase boson-boson interaction is frequently used here and in the literature, in the composite boson system (within the crossover 
scenario), these interactions should always be viewed as indirect and associated with the underlying fermionic interactions. 

In summary, the T — description of the composite boson problem (in the strong coupling limit) represents a cross between 
ideal and non-ideal Bose-Einstein physics. |-Whereas, an incomplete condensation at zero temperature is required in order to 
obtain superfluidity in a "true" Bose systen£j, in the crossover scenario, the condensation is always complete (at T — 0), but, 
nevertheless the phase mode and superconductivity exist. 

B. Physical Picture at Finite T < Tc: Present Work 

In this section we discuss our general results and physical picture for the finite T crossover problem, at an intuitive level. These 
results are then re-derived in more detail in Sec. Ill using a microscopic T-matrix approach. As the temperature is increased above 
T = 0, and the coupling g becomes sufficiently strong, two important effects ensueEJ. (i) The excitation gap is no longer the 
same as (the amplitude of) the superconducting order parameter, (ii) Incoherent pair excitations with non-zero center of mass 
momentum can be thermally excited. Indeed, points (i) and (ii) are inter-related. The first of these is anticipated as one leaves 
the BCS regime, where even above Tc, there is expected to be a normal state gap for fermionic excitatimis; this gap is associated 
with meta-stable (intermediate gjB or stable (large g) pairs. It is natural to presume, as we have foundjlj that these effects persist 
below Tc- Indeed, if there is an excitation gap at Tc, in any second order superconducting transition, this gap must necessarily 
be different from the superconducting order parameter, at least, at and slightly below Tc- 

The second of these two points was noted in Ref. nS. These pair excitations are related to the "particle" excitation branch 
of the interacting Bose sy.stem, which (although similar to collective phase mode branch at small wave-vectors), represents a 
distinct dynamical modeE3. In the composite boson problem the pair excitations and related pair propagator, called T{ifl, q), 
play an essential role. In the limit of non-zero Q = {iil, q), we refer to this propagator as Tpg{Q), where the subscript pg derives 
from "pseudogap". For small Vt and q, it may be approximated by 

Tpg{ifl, q) « Z/{in -flq + iTq + Upatr) , (5) 

where Z is the usual renormalization factor, ifi is a bosonic Matsubara frequency, fiq is the dispersion of the finite momentum 
pair excitations (with ilq=o = 0), F^^ is the pair excitation lifetime, and iJ,pair is the associated effective pair "chemical 
potential" with its value set by the well known (ideal gas) BEC condition 

^ipa^r^O, for T < Tc . (6) 

It is this function Tpg, which will be shown in Sec. Ill, to determine the fermionic excitation gap A(r) and chemical potential 
/i. Moreover, the fermionic excitation gap is related, in turn, to the condition that ^ipair = at and below Tc- Quite generally, 
the presence of incoherent pair excitations blocks the available states for fermions and, thereby, affects the fermion excitation 
gap. These effects, which occur both above and below Tc, are directly related to particle-hole asymmetry which appears as the 
system crosses out of the BCS regime. 



It should be stressed that within the BCS-BEC crossover picture, the form of T is highly circumscribed so as to produce the 
ground state equations, Eqs. ^ and (Q). Through these equations the fermionic degrees of freedom play an important role, even 
at very strong coupling. This constraint can then be used to deduce the function ilq. We find that the pairing approximation 
(discussed in Sec. Ill below) produces this T = description of the fermionic degrees of freedom, and that, as a consequence, 
for sufficiently small q, the pair excitation dispersion is given by 

17q - qV2Mp,i, (7) 

where the pair mass Afpaii is dependent on g, temperature, density, lattice structure and other materials properties. This dispersion 
is found to be a consequence of other T-matrix schemes, as well, and takes a similar form above Tc, although there the pairs are 
not as long lived. Once the pair propagator is characterized, Tc can be obtained either when approached from above or below. 

We may now quantify the first point (i) listed above. The deviation between the excitation gap and the order parameter is 
related to the number of thermally excited finite momentum pair excitations. We define the difference between the excitation 
gap A and order parameter A^c as 

Alg = A2 - aL , (8) 

where throughout this paper. Age is taken to be real. The number of (incoherent) pair excitations is given in terms of the pair 
propagator Tpg 
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(9) 



where h{Vt) is the Bose function and Tpg (fJ, q) is the analytically continued {iVt -^ il + iO+) pair propagator. In the regime of ia,- 
termediate coupling, where pseudogap effects are apparent, these bosonic excitations act in concert with fermionic excitationsM. 
This is what one expects if there is to be a smooth interpolation between the BCS and BEC limits. 

The behavior of Apg, A and A^c is schematically plotted in Fig. 1 in the three different regimes: weak (BCS), intermediate, 
and strong coupling (nearly BEC) regimes. Below T^ these plots are based on detailedjiumerical calculationslllj, whereas above 
Tc, where the computations are more difficultEJ, on a simple extrapolation procedures r, 

It may be noted that Fig. 1 has a direct analogue in the (Coulomb modulated) phase fluctuation scenarioH. The three panels 
(from top to bottom) would then correspond to progressively decreasing the size of the plasma frequency Up. In this phase 
fluctuation scenario the tuning parameter is LOp, whereas in the crossover scenario the parameter g sets the scale for the size of 
Ap . It should, thus, be clear that Coulomb-modulated phase fluctuations are not the only way to create an excitation gap which 
appears above the transition temperature. 

It should, finally, be stressed that these incoherent, finite momentum pair excitations, which enter into Apg via Eq. (^, are 
irrelevant in the BCS limit, in accord with Eq. (pi) and Fig. 1(a). In that limit the "quasi-particle" assumption implicit in Eq. (ph is 
invalid and because of both the lack of particle-hole asymmetry and the large damping Fq, the pair excitation spectrum overlaps 
the partiijle-particle continuum. One can quantify the reliability of this key approximation. In Fig. 2, using our numerical 
schemdiil, we plot the value A of the wave-vector |q| at which the pair dispersion intersects the continuum states, as a function 
of coupling g. This corresponds to a measure of the Landau damping of the pair propagator. The shaded region indicates where 
the incoherent finite momentum pairs represent ill-defined excitations. Outside this shaded region, the assumptions implicit in 
Eq. (Bl) should be valid. Related calculations show that Apg{Tc) / Asd^) is arbitrarily small in the weak coupling limit, so that 
even if we apply Eq. (||) directly in this limit, pseudogap effects are negligible in the BCS regime. 

C. The Strong Coupling Limit: "True" versus Composite Bosons 

A first important distinction can be found, between true and composite bosons, at the level of the Leggett ground state. 
As noted earlier, already at T = there appears to be a_jnix of quasi-ideal and interacting Bose gas character to the strong 
couphng limit. |-Xhje gap equation [Eq. (|])] is associatedEJ with "noninteracting diatomic molecules", whereas, the collective 
mode spectrumEjlD reflects an effective boson-boson interaction which relates to the Pauli statistics of the constituent fermions. 
This is revealed most clearly in jellium models where the AB sound velocity remains finite at infinite g, with an asymptote 
associated with the residual interactions. Some insight into the origin of these boson-boson effects can be found in Refs. |3^ and 
p5| . Indeed, one should be cautious in the use of the phrase effective boson-boson repulsion which, for large g, derives entirely 
from Fermi statistical effects. Unlike in the interacting Bose systerr£3, where boson-boson interactions need to be separately 
included in the boson propagator, here the physics associated with the Pauh principle is already accounted for and should not be 
fed back again to renormalize the dispersion of the composite bosons. 



A second important difference arises from the fact that this superconducting ground state corresponds to one in which there is 
full condensation so that, as in the BCS phase, the condensate fraction no = n. By contrast, in a Rose superfluid there is always 
a depletion of the condensate at T = 0, caused by the existence of the boson-boson repulsion. 

As a final important difference, we note that the behavior of the pair propagator T, which must necessarily be consistent with 
Eqs. (Bj) and &), is highly circumscribed and rather different from what one might deduce based on the Bogoliubov model for 
a Bose liquidH. The fermionic degrees of freedom can never be fully "integrated out". The fermionic excitation gap A and the 
pair chemical potential /ipair = are, moreover, closely inter-related via, e.g., Eq. ( p^ below. These effects have no natural 
counterpart in the Bose liquid (where the fermionic excitation gap is of no consequeac^ JU is, thus, not surprising that in Tc 
calculations, which are associated with a divergence in T, a variety of different groupstj lScJ find that the strong coupling limit 
is characterized by a quasi-ideal BEC result 



Tr^ 



3.31 



M, 
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air ^ ^ ' 
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(10) 



where A/, 



pair is the same pair mass as that which appears in Oq = c^ /Mpair of Eq. (^. Thus, from the perspective of Tc 



and of the temperature dependent gap pa 
have reached rather similar conclusionsJlS' 



ns, the bosons are "free", except for the renormalized mass. Related theories 
^' These conclusions owe their origin to the underlying mean field structure of BCS 



theory, which is the starting point of the crossover scenario, 
general T-matrix approaches, as is discussed below. 



At a more microscopic level, they appear to be associated with 



III. SELF-CONSISTENT T-MATRIX APPROXIMATIONS 



A. General Results for the Superconducting Phase 



In this section we present the self consistency conditions and related gap equations associated with the superconductingstate 
within the broad class of T-matrix-based crossover theories. An important goal of this discussion is to show that Eqs. (ot-(H) 
are rather general consequences of these schemes when applied below Tc- In these non-perturbative approaches the two-particle 
Green's function is expressed in terms of a T-matrix (or pair propagator), T{Q), which is determined self consistently in terms 
of the single particle Green's function. One solves, in effect, three coupled equations for the self energy, T-matrix and chemical 
potential (via the number equation): 



^{K) - G-\K) - G~\K) = Y. nQ) G{Q - K) 



(11a) 



9 = [1+3X(Q)]T(Q) 



(lib) 



2Y^G{K) 



K 



(lie) 



Here and in what follows we use four-vector notation: Q = (ifl, q), ^g = ^Sin Sq' ^'^e. The choice of the functions G 
and pair susceptibility x varies from one approximation to another Here we address two-different schemes which introduce 
self consistency at a level beyond the lowest order T-matrix approximation used by NSR,ta so that our discussion focuses on 
T-matrix schemes where dressed GreeuI^^£ipctions enter into the self consistency requirements, (a) The first approach, known 
as the FLEX (or 'GG') approximationLJ'Ej'O takes G — G with x = x^^^'^ where we define 



FLEX 



{P) =Y,G{K)G{P - K) 



(12) 



K 



with the corresponding self energy 



E^^^-^(i^) = Yr{P)G{P - K) 



(13) 



In addition, we study the (b) pairing approximationEBM (or 'GGo 'scheme), which sets G = Go and x. = x^°'^^ with 

;^P--(P)=^G(i^)G,(F-i^) 



(14) 



with the corresponding self energy 

EP"''-(if)^^r(P)Go(P-if) (15) 

p 

We know of no other literature on the FLEX scheme within the superconducting state (appropriate to a fully three dimensional 
system). By contrast, abovCpTc there is considerable literature on the behavior of T and Tc, which, in the FLEX scheme, have 
been calculated numericallyL3 and analyticallyt^ by solution of Eqs. i\\A\). As the temperature is lowered, this T-matrix develops 
a maximum around Q = and the transition temperature Tc to the broken symmetry phase is signaled by a pole given by the 
condition 

gT-\Q = 0;Tc) = l+gx{Q = 0;Tc,A)=0. (16) 

The counterpart of this analysis is considerably more complicated below Tc- The procedure, which we summarize in this 
section, is an approximation which is chosen to satisfy the following four criteria: (i) It leads to the same transition temperature 
when approached from below as from above, (ii) It leads directly to a physical interpretation of the BCS-BEC crossover scheme 
as discussed in the previous section and illustrated in Fig. 1. An important, and third criterion which we view as an additional 
check on the approximations used is that (iii) one should recover the BCS scheme in weak coupling for all T < Tc and, finally, 
(iv) we should recover the Leggett ground state at T = 0. 

An important premise underlying these crossover schemes is that the particle-particle dominates the particle-hole channel. 
Without such an assumption one would be forced to evaluate a full 2x2 Nambu matrix Green's function along with a 4 x 4 
T-matrix. A complete self-consistent solution of the resulting set of equations is prohibitively difficult. With the assumption that 
the particle-particle channel is dominant, one argues that formally the set of Eqs. ( |ll| ) remains valid below Tc, provided that the 
T-matrix acquire a singular delta-function component, which describes the (Q = 0) Cooper pair condensate in equilibrium, so 
that 

T{Q)=Zc{Q)+Tpg{Q), (17a) 

where 



fsiQ) , and r,,iQ) ^ ^^^^ . 



%c{Q) = -^S{Q) , and TpgiQ) = ^^-j^ . (17b) 



This form for T guarantees that the regular part Tpg of the T-matrix remains finite for any non-zero Q = {ifl, q). Inserting 
Eqs. (HTh into Eq. (lib) and treating separately the delta function contribution and regular terms, one arrives at the following 



quite general gap equation 

l + 5X(Q = 0;r,A) = 0, T<Tc. (18) 

In our self consistent scheme, the critical temperature Tc can be obtained from Eq. (jlj) by setting Asc, (which is implicitly 
contained in A), to zero. This result coincides precisely with the Thouless criterion of Eq. (|l6|), which is obtained by approaching 
the transition from above. As a result of the form of Eq. (|I8|), the regular part of the T-matrix Tpg (Q) diverges as Q -^ 0; it can 
be written in the form 

Tpg{in,q) = z/itn-n^ + irc), t<Tc, (19) 

in accord with Eq. (|^. Here Fq -^ as q ^ 0. Equation (18) also leads directly to a microscopic derivation of Eq. (g). The self 



energy of Eq. (11a) may be decomposed into two terms: 

S(X) = ^sc{K) + Epg{K), (20) 

where the term associated with the condensate contribution, called Tgc is 

^sciK) - -^i.Gi-K) . (21) 

In evaluating the pseudogap contribution to S, as a consequence of Eq. (Il8|), the main contribution to the Q sum comes from the 
small Q region, so that the integral may be approximated by 

^pgiK) « Gi-K) Y^TpgiQ) = -AlgG{-K) . (22) 

Q 



In this way the total self energy is 



i:{K) w -A'^G{-K) , A 



A2 +A2 

sc ' pg 



(23) 



Here the total excitation gap is to be associated with the parameter A, and we see that Eqs. (|2|) and (^3|) naturally lead back to 
the definition of Apg which appears in Eq. (Q). 

There is another way of relating A^ to the fluctuations of the pairing field, which leads to a more precise interpretation of 

Eq. (|^). We write the pairing field Aq(t) = \g\ X]k'^-k+q/2i(^)ck+q/2T(^)' about its mean field value (|A|) = |(Aq=o(0))| = 
A,,. We define (|A|2) = Y.q{^Q^q) = 9^Y.k,k'Y.qC2{K,K'-Q), where C2{K,K'-Q) is the proper two-particle 
correlation function, which in general is not factorizable in the variables K and K'. After some algebra one arrives at 



(|Ap)=-^T(Q)[5x(Q)f 
Q 



Now, inserting in Eq. ( |24[ ) the expression for the T-matrix given by Eqs. ( |17[ ), we obtain 

l2 



\A\')=Al-J2T,,{Q)[gx{Q)r 
Q 



(24) 



(25) 



Since Tpg is highly peaked about Q ~ 0, the pair susceptibility xiQ) on the right hand side in Eq. (^SJ) can be approximated by 
its Q = value. We again make use of the gap equation ([iq) to write y^^ 7pg(Q)[.9x(Q)]^ ~ ^^o'^pgiQ) — ^pg- Thus, 
quite generally 



A2 

pg 



Q 



(26) 



The above discussion is expected to apply to both the FLEX scheme and the pairing approximation. Moreover, on the basis of 
the behavior above Tc, there is no a priori reason to select one approach over the other However, the latter seems to be preferred 
if one imposes the third and fourth criteria discussed above. 

Indeed, the superconducting state is associated with the self energy of Eq. (|2l|). It can be seen that E^c coincides with the BCS 
self energy if we adopt the pairing approximation so that G = Go- With this result the BCS gap equation follows from Eq. (Uq): 

1 + g ^ G{K)Go{-K) = 1 + g ^ ,^^^^, ^ . (27) 



Q 



K 



K 



The results of the pairing approximation, at all g, can be summarized simply. This scheme leads to the following equations 
for the three unknowns A , /^, and Apg 



9 



-^ + E 



1 - 2fiEk) 
2Ek 







(28a) 



where Ei^ = y'e^ + A^. This equation must be solved self consistently with 



2E 



I + ^ fiE^) 



where t;2 ^ i(l - et/^k). 

Finally, the decomposition of A into Age and Apg requires the solution of a third equation, namely Eq. 
here for completeness 



/oo HO 
— b{n)im%g{n,ci) 



(28b) 

which we repeat 
(28c) 



An additional che ck o n the validity of the approximation scheme relates to criterion (iv), that is, the nature of the ground state. 
It follows from Eq. (28c) that, as a result of the Rose function &(57), quite generally 



lim A 



pg — L) , 



(29) 



as is consistent with a physical picture in which Apg is associated with classical fluctuations. When this limit is used within the 
pairing approximation, the ground state which results from Eqs. (28a) and (28b) is the same as that proposed by Leggett [see 
Eqs. (|)and(H)]. 



IV. ELECTROMAGNETIC RESPONSE AND COLLECTIVE MODES OF A SUPERCONDUCTOR: BEYOND BCS THEORY 

The purpose of this section is to study the gauge invariant (linear) response of a superconductor to an external electromagnetic 
(EM) field, and obtain the associated collective mode spectrum. Our discussion is generally relevant to complex situations 
such as those appropriate to the BCS-BEC crossover scenario. An important ingredient of this discussion is establishing the 
role of particle-hole asymmetry. It should be noted that there.|aie fairly extensive discussions in the literature on theJaehavior 
of collective modes within the T — crossover scenarioEJ'LJt^. Here we review a slightly different formulationEZrE3 which 
introduces a matrix extension of the Kubo formalism of the normal state. We find that this approach is more directly amenable 
to extension to finite T, where the pair fluctuation diagrams need to be incorporated. 

The definition of the collective modes of a superconductor must be made with some precision. We refer to the underlying 
Goldstone boson of the charged or uncharged superconductor as the "AB mode", after AndersonEj and BogoliubovcJ. According 
to our specific definition, this AB mode appears as a pole structure in the gauge invariant formulatiops of the electrodynamic 
response functions, for example, in the density-density correlation function. Early work by PrangeEJ referred to this as the 
"ghost mode" of the neutral system, since this term is not directly affected by the long range Coulomb interaction. By contrast, 
the normal modes of the charged or uncharged superconductor, which we shall call the "collective modes", involve a coupling 
between the density, phase, and for the BCS-BEC case, amplitude degrees of freedom. For these, one needs to incorporate a 
many body theoretic treatment of the particle-hole channel as well. In crossover theories this channel is not as well characterized 
as is the particle-particle channel. 

A. Gauge Invariant EM Response Kernel 

In the presence of a weak externally apphed EM field, with four-vector potential A'' ~ (0, A), the four-current density 

J'^ = (p, J) is given by 

J^{Q) = K'^''{Q)A,{Q), (30) 

where, Q = q'^ — {u!,q) is a four-momentum, and K'^'^ is the EM response kernel which can be written as 

K'^-'iQ) = K^'iQ) + SK^^^iQ) . (31) 

Here 

is the usual Kubo expression for the electromagnetic response. We define the current-current correlation function P^'^{t, q) = 
— i6'(r)([j'^(T,q), j'^(0, — q)]). In the above equation, g^'' is the contravariant diagonal metric tensor, with diagonal elements 
(1, —1, —1, —1), and n, e and m are the particle density, charge and mass, respectively. 

The presence of 5K^^'^ in Eq. (|T]) is due to the perturbation of the superconducting order parameter by the EM field, i.e., to 
the excitation of the collective modes of lS.sc- This term is required to satisfy charge conservation g^ J^ — 0, which requires that 

q^.K^'^iQ) - . (33a) 

Moreover, gauge invariance yields 

K^"{Q)qu = , (33b) 

Note that, since K^^^{—Q) — K'^^{Q), the two constraints Eqs. ( P3| ) are in fact equivalent. 

The incorporation of gauge invariance into-a general microscopic theory may be implemented in several ways. Here we do so 
via a general matrix linear response approachE3 in which the perturbation of the condensate is included as additional contributions 
Ai + i/S.2 to the applied external field. These contributions are self consistently obtained (by using the gap equation) and then 
eliminated from the final expression for K^^'^ . We now implement this procedure. Let 771 2 denote the change in the expectation 
value of the pairing field ^1.2 corresponding to Ai_2. For the case of an s-wave pairing interaction g < 0, the self-consistency 
condition Ai 2 = ~ff'7i,2 leads to the following equations 

jp = K^"'A^ = K^^A^ + ff'^Ai + i^^^^a , (34a) 



m = -— = R'"A, + QnAi+Qi2A2, (34b) 

g 

V2 = -— = R^'^A, + Q21A1 + Q22A2 , (34c) 
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where i?^*(r,q) = -i6l(T)([j''(T, q),7?,(0, -q)]), with ^ = 0, . . . , 3, and i = 1,2; and Qij(r,q) = 
-i9{T){MT, q), %(0, -q)]), with ^, J == 1, 2. 

Thus far, the important quantities -K^^", i?^* and Qij are unknowns which contain the details of the appropriate microscopic 
model. We shall return to these later in Sec. ^ The last two of Eqs. ( p4| ) can be used to express Ai 2 in terms of A^, 

Q22R " ~ Q12R " . ,,- . 

Ai = --= — = A^ , (35a) 

Q11Q22 — Q12Q21 

QiiR " - Q21R " , ,T,, , 

A2 = — = — = A^ , (35b) 

Q11Q22 — Q12Q21 



where Qu = 1/g + Qu, with i — 1,2. Finally, inserting Eqs. (35) into Eq. (34a) one obtains 



with 

^^^. ^ QiiR^^R^" + Q22R'''R"' ~ Qi2R^^R^'' - Q2iR''^R"' ^3^^^ 

Q11Q22 — Q12Q21 

As can be seen from the above rather complicated equation, the electromagnetic response of a superconductor involves many 
different components of the generalized polarizability. Moreover, in the form of Eqs. ( [36b| ) it is not evident that the results are 
gauge invariant. In order to demonstrate gauge invariance and reduce the number of component polarizabilities, we first rewrite 
K'^" in a way which incorporates the effects of the amplitude contributions via a renormalization of the relevant generalized 
polarizabilities, i.e., 

j^t,^ = K'^" + 6K""" , (37a) 

where 

^ = -^0 ^ ' (37b) 

yii 



and 



R"^' = R^^-^R'^\ Q'^^^Q,,-92^, (37c) 



In this way we obtain a simpler expression for 6K 



11 Qii 

,fJ,U 



p/^2 T>l2i' 

SK'"" = - , . (38) 



' 22 



We now consider a particular (a priori unknown) gauge A'^ in which the current density can be expressed as J^ = K"^f'^A'^. 
The gauge transformationcd which connects the four-potential A^ in an arbitrary gauge with A' i.e.. A' = Afj, + ixqfj., must 
satisfy 

J'^ = K'^-'A, = K'^" {A, + ixq.) . (39) 



Now invoking charge conservation, one obtains 






and, therefore, 

J^^,u ^ ^,^'> _ A 1\. L . (41) 

The above equation satisfies two important requirements: it is manifestly gauge invariant and, moreover, it has been reduced to 
a form that depends principally on the four-current-current correlation functions. (The word "principally" appears because in 
the absence of particle-hole symmetry, there are effects associated with the order parameter amplitude contributio ns w hich enter 
via Eq. ( pTI ) and add to the complexity of the calculations). Equation (^Tj) should be directly compared with Eq. (36b). In order 



for the formulations to be consistent and to explicitly keep track of the conservation laws (p3h, the following identities must be 
satisfied: 

(gM<'^)Q'22 = ('ZM^'^')^''^ (42a) 

These identities may be viewed as "Ward identities" for the superconducting two particle correlation functionsc3. Any theory 
which adds additional self energy contributions to the BCS scheme must obey these important equations. We shall return to this 
issue in Sec. 0. 

B. The Goldstone Boson or AB Mode 

The EM response kernel [cf. Eqs. (^7[)-(p4|)] of a superconductor contains a pole structure which is related to the underlying 
Goldstone boson of-the system. Unlike the phase mode component of the collective mode spectrum, this AB mode is independent 
of Coulomb effectsEJ. The dispersion of this amplitude renormalized AB mode is given by 

q^.K^'^, = . (43) 

For an isotropic system K'^ — K'^^ 5ap, and Eq. ( ^3| ) can be rewritten as 

cj'K'''„" + eKf-2u;q^K''o" = 0, (44) 

with a = 1,2,3, and in the last term on the LHS of Eq. (Q) a summation over repeated Greek indices is assumed. It might 
seem surprising that from an analysis which incorporates a complicated matrix linear response approach, the dispersion of the 
AB mode ultimately involves only the amplitude renormalized four-current correlation functions, namely the density-density, 
current-current and density-current correlation functions. This result is, nevertheless, a consequence of gauge invariance. 

At zero temperature K'q" vanishes, and the sound-like AB mode has the usual linear dispersion uj = ujq ^ c|q| with the 
"sound velocity" given by 

c'^Kf/K'l\ (45) 

The equations in this section represent an important starting point for our numerical analysis. 

C. General Collective Modes 

We may interpret the AB mode as a special type of collective mode which is associated with Aiy = Oin Eqs. (p5|). This mode 
corresponds to free oscillations of Ai 2 with a dispersion lo ^ cq given by the solution to the equation 

detlQy I = Q11Q22 - Q12Q21 = . (46) 



More generally, according to Eq. (34a) the collective modes of the order parameter induce density and current oscillations. In 
the same way as the pairing field couples to the mean field order parameter in the particle-particle channel, the density operator 
p{Q) couples to the mean field 6(t){Q) — V{Q)6p{Q), where V{Q) is an effective particle-hole interaction which may derive 
from the pairing channel or, in a charged superconductor, from the Coulomb interaction. Here Sp = (p) — po is the expectation 
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value of the charge density operator with respect to its uniform, equilibrium value po- Within our self-consistent linear response 
theory the field S(f) must be treated on an equal footing with Ai 2, and formally can be incorporated into the linear response 



'MO s 



of the system by adding an extra term Kq 5(j) to the right hand side of Eq. (34a). The other two Eqs. ( |34| ) should be treated 
similarly. Note that, quite generally, the effect of the "external field" (50 amounts to replacing the scalar potential ^° = </) by 
A'^ = 4) = 4) + d(j). In this way one arrives at the following set of three linear, homogeneous equations for the unknowns S^, Ai, 
and A2 

= i?i",50 + QnAi + Q12A2 , (47a) 

= R^"d(f> + Q21A1 + Q22A2 , (47b) 

Sp^^-^ = K^o^'S^ + i?"iAi + R"'A, . (47c) 

The dispersion of the collective modes of the system is given by the condition that the above Eqs. (B7h have a nontrivial solution 

Qii + 1/5 Q12 R^° 

Q21 Q22 + 1/5 R^° 

^01 ^02 ^00 _ ^/y 



= . (48) 



In the case of particle-hole symmetry Q12 ~ Q21 = R^'^ = R'^^ = and, the amplitude mode decouples from the phase and 
density modes; the latter two are, however, in general coupled. 

V. EFFECT OF PAIR FLUCTUATIONS ON THE ELECTROMAGNETIC RESPONSE: SOME EXAMPLES 

Once dressed Green's functions G enter into the calculational schemes, the collective mode polarizabilities (e.g., Q22) and the 
EM response tensor K^" must necessarily include vertex corrections dictated by the form of the self energy E, which depends 
on the T-matrix T, which, in turn depends on the form of x- These vertex corrections are associated with gauge invariance 
and with the constraints which are summarized in Eqs. (Q). It can be seen that these constraints are even more complicated 
than the Ward identities of the normal state. Indeed, it is relatively straightforward to introduce collective mode effects into 
the electromagnetic response in a completely geneijal fashion which is required by gauge invariance. This issue was discussed 
in Sec. IV as well as extensively in the literatureElO. The difficulty is in the implementation. In this section we begin with a 



discussion of the T = behavior where the incoherent pair excitation contributions to the self energy corrections and vertex 
functions vanish. In this section, we shall keep the symmetry factor (^k explicitly. 

A. r = Behavior of the AB Mode and Pair Susceptibility 

It is quite useful to first address the zero temperature results since there it is relatively simple to compare the associated 
polarizabilities of the AB mode with that of the pair susceptibility x- In the presence of particle hole symmetry this collective 



mode polarizability can be associated with Q22 which was first defined in Eq. (34c). In the more general case (which applies 
away from the BCS limit) Q22 must be replaced by a combination of phase and amplitude terms so that it is given by Q22 = 
Q22 — Qi2Q2i/Qii- 

We may readily evaluate these contributions in the ground state, where Asc = A. The polarizability Q22 is given by 

Q22(g) ^^Y. [Gi~P)G{P -Q) + G{P)G{Q -P)+ F^{P)F\P -Q) + F{P)F{P - Q)] ^l^^^^, (49) 

^ p 

where 

G{K) = Go{K)/[l + Al^iplGohK)GoiK)] , (50) 

and 

F{P) = AsciPpGiP)Goi-P) . (51) 

Now, it can be seen that the pair susceptibility x in the pairing approximation satisfies 

J2[G{-P)GiP) + F{P)F{P)]vl = Y. G{P)Go{-P)vl = X(0) (52) 

p p 
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and, moreover, Qi2(0) = Q2i(0) = so that 

- + g22(0) = i[l + <?x(0)] = 0. (53) 

9 9 

In this way, the AB mode propagator is soft under the same conditions which yield a soft pair excitation propagator Tpg = 
g/{l + 5x)' ^nd these conditions correspond to the gap equation Eq. O). Moreover, it can be seen that Q^iiff) = Q'2i{ Q) so 
that, upon expanding around Q = 0, one has (322 (Q) = "22^^-/922? , Qi2(Q) = -Q2i(Q) = -i^a^^, andl/5 + (3ii((3) = 
l/.g + ail, where 



/322 - ;7 E ?? 






^k(V6k)' - 7(v4) ■ (v^l) 



Etk 2 
k k 

i:4^^ (54) 



an = 

where d denotes the dimensionality of the system. Thus, one obtains 

c = ^2 . (55) 

«22 + 17^:^7 

At weak coupling in 3D, where one has particle-hole symmetry, ai2 = 0, the amplitude and the phase modes decouple. This 
leads to the well known result c ^^v^l\fz. More generally, for arbitrary coupling strength g, these equations yield results 
equivalent to those in the literaturecjOEJ, as well as those derived from the formalism of Sec. IV.B. Finally, it should be noted 
that, since both Eqs. (Wlh and ( p^ ) have the same poles, the condition Q^iiff) = yields the same AB mode dispersion as that 
determined from Eq. ([43|). This is a consequence of gauge invariance. 

B. AB Mode at Finite Temperatures 

We now turn to finite temperatures where there is essentially no prior work on the collective mode behavior in the crossover 
scenario. At the level of BCS theory (and in the Leggett ground state) the extended "Ward identities" of Eqs. ( |44 ) can be 
explicitly shown to be satisfied. Presumably they are also obeyed in the presence of impurities, as, for example, in the scheme of 
Ref. |3^. However, in general, it is difficult to go beyond these simple cases in computing all components of the matrix response 
function. Fortunately, the calculation of the AB mode is somewhat simpler. It reduces to a solution of Eq. (pi), which, in the 
presence of particle-hole symmetry, involves a computation of only the electromagnetic response kernel: the density-density, 
density-current and current-current correlation functions. 

It is the goal of this section to compute these three response functions within the "pairing approximation" to the T-matrix. 
Our work is based on the normal state approach of PattonEZi and the associated diagrams are shown in Fig. 3. Because full 
Green's functions G appear in place of Go (as indicated by the heavy lines) these, diagrams are related to but different from 
their counterparts studied by Aslanwp¥ and Larkin and by Maki and ThompsonE3. This diagram scheme forms the basis for 
calculations published by our groupll3lla of the penetration depth within the BCS-BEC crossover scheme. 

Here we make one additional assumption. We treat the amplitude renormalizations which appear in Eqs. (pTI), only approxi- 
mately since these contributions introduce a variety of additional correlation functions, which must be calculated in a consistent 
fashion, so as to satisfy Eqs. (p^. Because the amplitude mode is gapped, at least at low T, we can approximate these amplitude 
renormalizations by their T = counterparts, which are much simpler to deduce. 

The three electromagnetic correlation functions reduce to a calculation of P^^ , which can be written as 

P>'''{Q) = 2 ^ A^(if, i^ - Q)G{K)G{K - Q)k''{K, K - Q) , (56) 

K 

where X{K, K - Q) = {I, ^jf^) and A{K, K - Q) ^ X{K, K - Q) + 6Asc{K, K - Q) + 6Apg{K, K ~ Q) aie the bare 
and full vertices, respectively. 
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To evaluate the vertex function A'' we decompose it into a pseudogap contribution Ap^ and a superconducting contribution 
Age- (The latter can be regarded as the Gor'kov "i^" function contribution, although we do not use that notation here). The 
pseudogap contribution comes from a sum of Maki-Thompson (MT) and Aslamazov-Larkin (ALi 2) type of diagrams [see 
Fig. 3(b)]. Since these vertex corrections can be obtained from a proper vertex insertion to the self energy, it follows that 
there is a cancellation between these various terms which simplifies the algebra. This cancellation is shown in more detail in 
Appendix 0. Following the analysis in this Appendix, the sum of both pg and sc contributions is given by 



SA^'{K,K-Q)^~{A 



pg' 



AL)^k<^k-qGo(-/OGo(0 - K)X''iQ - K, -K) 
AL Go (-if )-"'''''"''/' 



d^l- 



pg 



dka 



(57) 



where use has been made of the fact that Tpg (Q) is highly peaked at Q = 0, and that A'i = 



^Q ■'pgyQj- 



The AB mode dispersion involves the sum of three terms which enter into Eqs. (M3fi and (M4b. We next substitute Eqs. (BTh 
into Eq. (56|. After performing the Matsubara frequency summation, and analytically continuing iil ^ il + iO+, we obtain for 
small J7 and q 



q^K^^Qu = q • ( - + P ) • q - 2r!q • Po + n'Pf 



00 






-2^^E 



1 - 2/(^k) 



2Ek 



I'{E^) 



^k(Vekf-4(V6^)-(V(p£) 



c'Pk 



El 



' - '^(^-^^ f fiE.) - fiE,)^- (^ ■ ^'^^' - ^'(^ ■ ^^^^' 



2^k 



172 - (q . VEk)^ 



,2\2 



-^/iE.) 



(q-v4)(q-vyk) + A^(q- y^l) 

f72 - (q • VEk)'^ 



(58) 



where f{E) is the Fermi function. Because Eq. ( pq ) is ill-behaved for long wavelengths and low frequencies, in order to calculate 
the AB mode velocity one needs to take the appropriate limit ft — cq ^ 0. By contrast, the calculation of the London penetration 
depth first requires to set 57 — (static limit), and then q ^- 0. The superfluid density rig can be calculated from the coefficient 
of the q2 term in Eq. (pSb, (see, also Eq. (Bl ) for Q = 0). Finally, the AB mode "sound" velocity c — fl/q, in the absence of 



the amplitude renormalization, can be obtained by solving q^Kl^ g^ = 0. 

In the absence of the pseudogap (i.e., when A^^ = lA) the last term inside {. . .} in Eq. ( ^8] ) drops out, and the resulting 
analytical expression reduces to the standard BCS resulto, which at T = has the relatively simple form 



^2 a2 

k 



d 4- El 



^k(V6kf-7(v4) 



(Vy^^) 



^'E 
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cVl 



El 



(59) 



At finite T, the AB mode becomes damped, and the real and imaginary parts of the sound velocity have to be calculated 
numerically. Although, the algebra is somewhat complicated, it can be shown that the AB mode satisfies c — > as T -^ r^, as 
expected. 

VI. NUMERICAL RESULTS: ZERO AND FINITE TEMPERATURES 

In this section we summarize numerical results obtained for the AB mode velocity c associated with the electromagnetic 
response kernel, as obtained by solving Eqs. ( |43| ) and (^^. We also briefly discuss the behavior for the T — phase mode 
velocity v^ which results from the coupling to density..|fluctuations, as well [see Eq. (Sq)]. The former, which has physical 
implications for the behavior of the dielectric constantMEa, is the more straightforward to compute, because it does not require 
any new approximations associated with the effective interactions V in the particle-hole channel. The analysis of this section 
provides information about the nature of the "quasi-ideal" Bose gas limit, which we address via plots of the infinite g asymptote 
of the AB mode, called Coo- It also helps to clarify how pair fluctuations contribute, at finite temperatures, to the collective mode 
dispersion. Our T = calculations are based on the Leggett ground state which corresponds to that of the pairing approximation 
as well. At finite T, we numerically evaluate the AB sound dispersion from Eq. (pg), obtained within the framework of the 
pairing approximation. 

In Fig. 4 we plot the zero temperature value of c as a function of the dimensionless coupling strength g/gc, where gc — 
—4:Tr/mko is the critical value of the coupling above which bound pairs are formed in vacuum. Here we consider a 3D jellium 
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model at three different electron densities, which are parameterized via ko/kp. The most interesting feature of these and 
related curves is shown in the inset where we plot the asymptotic limit for each value of density or, equivalently, ko- This 
numerically obtained asymptote reflects the effective residual boson-boson interactions in the "quasi-ideal" Bose gas limit, 
and is close to the value calculated in Ref. ^ whose functional dependence is given by Coa/vp oc ^/kp/ko or, equivalently. 
Coo fx Y^n/fco- Interpreting the physics as if the system were a true interacting Bose system, one would obtain an effective 
interaction C/(0) « Sn"^ /mko, independent of g in the strong coupling limit. As expected, these inter-boson interactions come 
exclusively from the underlying fermion character of the system, and can be associated with the repulsion between the fermions 
due to the Pauli principle. All of this is seen most directlyta by noting that the behavior displayed in the inset can be interpreted 
in terms of the effective scattering length of the bosons as which is found to be twice that of the fermions ap in the strong 
coupling limit. Effects associated with the coupling constant g are, thus, entirely incorporated into making bosons out of a 
fermion pair, and are otherwise invisible. 

The same calculations are repeated in Fig. 5 for a tight binding lattice bandstructure at T = 0. Figure 5(a) plots the sound 
velocity for different densities n, as a function of the coupling constant; the behavior of the large g limit is shown in Fig. 5(b) 
as a function of density for a fixed g. Near half filling, where there is particle-hole symmetry, the amplitude contributions are 
irrelevant and the large g limit for c, which follows from Eq. (p9|), is c = \/2^, where t is the hopping integral. At low n the AB 
velocity varies as ^/n, which is consistent with the results shown above for jellium. In both cases the behavior again reflects the 
underlying fermionic character, since it is to be associated with a Pauli principle induced repulsion between bosons. Unlike in 
the jellium case, where c approaches a finite, asymptote as g increases, here c vanishes asymptotically due to the increase of the 
pair mass associated with lattice effectsJl3Ej For completeness, we also show, as an inset in Fig. 5(b), the behavior of v^, where 
we have used the RPA approximation to characterize the parameter V in the particle-hole channel. This approximation is in the 
spirit of previous work by Belkhir and RanderiaEj, although it cannot be readily motivated at sufficiently large g. 

Finally, in Fig. 6 we plot the temperature dependence of the AB mode velocity (both real and imaginary parts), for moderately 
strong coupling (solid lines) and the BCS limit (dashed lines). For the former, our curves stop somewhat below Tc, since close 
to the critical temperature, inaccuracies are introduced by our neglect of the temperature dependence of the amplitude mode 
contributions to c. It should be noted that the transition from a finite to zero value for the "sound" velocity appears to be rather 
abrupt in the vicinity of Tc, the stronger the coupling. This figure suggests that the AB mode velocity reflects the same transition 
temperature Tc as is computed via the excited pair propagator, or T-matrix. This represents an important self consistency check 
on the present formalism. 

VII. CONCLUSIONS 

This paper deals with the fairly complex issues of pair fluctuations, collective modes and gauge invariance in a BCS Bose- 
Einstein crossover scenario. A starting point for our approach is the Leggett ground state, which imposes rather strong constraints 
on the nature of the physics of fermions and composite bosons. The fermion degrees of freedom are always present through the 
self consistency conditions and can never be fully integrated out. Not only are these fermion pairs different from true bosons, but 
they represent a very special type of composite boson. They behave rather specifically, in correspondence with the underlying 
structure of the BCS state. E\pji at T = one_can see from previous work on the BCS-BEC crossover, that these constraints 
lead to a mix of ideal Bose gasEJ and non-idealo Bose liquid behavior. This mirrors some of the effects of BCS theory, in which 
the system undergoes a form of Bose condensation with a full condensate fraction rio/n — 1, at T = 0. Nevertheless, a BCS 
superconductor has a sound-like collective (phase) mode which is intimately associated with its superconductivity. Moreover, 
as is discussed in Sec. V.A, it is instructive to contrast the polarizability associated with the phase mode, called Q22{Q), with 
the pair susceptibility x{Q)- These two modes correspond to distinct dynamical branches, although both become soft at q = 
under the same conditions. The softness of the former is naturally associated with the Goldstone boson and the latter with a 
vanishing chemical potential for pairs: ^pair = 0. 

A central physics theme of this paper is reflected in Fig. 1. The crossover problem introduces an important new parameter 
Apg which characterizes the excited states of a non-BCS superconductor. As a result, there are three coupled equations to be 
satisfied for Apg, A and fi, at finite T, in contrast to zero temperature crossover theories where there are only two. This new 
parameter is a measure of the difference between the excitation gap and the superconducting order ppameter One can also 
arrive at a picture which is similar to Fig. 1 within the Coulomb-modulated phase fluctuation scenarioB. What is different here 
is the "tuning parameter", which corresponds to the coupling strength g in the crossover picture, and the plasma frequency ujp in 
the phase fluctuation scenario. 

One may summarize our results by asking the following series of questions, which our paper raises and answers. 

Is there any precedent for soft modes other than Goldstone bosons? 

Yes. In the present paper we find that a vanishing value for the pair chemical potential fipair below Tc leads to a soft mode 
corresponding to incoherent, finite momentum pair excitations. These excitations are analogous to the "particle" excitations in 
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the case of the neutral Bose hquid. Moreover, for the Bose hquid, this branch is also soft and always different (except at zero 
wavevector) from the Goldstone boson. However, in the true Bose system the two branches have the same slope at q = 0. 

Should the T-matrix T be renormalized so as to yield the sound mode dispersion at small wave-vectors, as in a Bose liquid? 
With this renormalization T would then be similar to its Bose liquid counterpart^ and in this way the pair excitation dispersion 
ilq at small wavevectors would be linear in q. We answer this question by noting that the fermion degrees of freedom strongly 
constrain T so that the composite boson system is different from the Bose liquid. As a consequence the proposed renormalization 
seems problematical for two reasons. Adding in this collective mode effect associated with the Bose liquid is equivalent to 
including boson-boson interactions. In the present composite boson case, these boson-boson interactions derive from fermionic 
degrees of freedom, i.e., the Pauli principle, which have already been accounted for in our calculations of T. It is not clear 
why these boson-boson interactions should then be included yet a second time. Secondly, once there is a renormalization of T, 
this will change the ground state gap equation and corresponding constraint on the fermionic chemical potential /i, associated 
with the changed fermionic self energy S. Moreover, this renormalization will, presumably, introduce an unphysical incomplete 
condensation in the ground state, like its counterpart in the Bose liquid. In summary, the composite boson system is considerably 
different from the true boson system, because the underlying fermionic constraints (through A, /i) can never be ignored. 

What about Coulomb renormalizations of the pair propagator? 

If Coulomb interactions were included, it would follow, if the analogy were appropriate, that the pair fluctuation mode would 
then be gapped. This would again compromise the self consistent conditions which must be satisfied from the fermionic per- 
spective (i.e., A and ^) in the Leggett ground state. Indeed, within the BCS formalism (as well as in the Leggett ground state) 
long range Coulomb interactions do not enter in an important way to change the gap equation structure, but rather they princi- 
pally affect the collective modes. It is for this reason that we argued earlier that Coulomb effects are presumed to be already 
included in the pairing interaction. It should be recalled that at large g, we have seen that essentially all signs of the two body (g 
dependent) fermion-fermion interaction are absent in the effective boson-boson interaction. 

Does this paper in any way change the way we think about BCS theory? Absolutely not. BCS theory appears as a special case 
in the weak coupling g, particle-hole symmetric limit of our more general approach. When we study pair excitations in this limit, 
we find they are greatly damped at all wavevectors q, and it makes no sense to talk about them. In this way Apg(Tc)/Asc(0) is 
vanishingly small in the BCS limit. 

To what extent are the results of this paper limited by the T matrix approximation ? The T-matrix approximation seems to 
be intimately connected to the physics of the crossover scenario. This scheme represents, in some sense, a truncation of the 
interactions at a pair-wise level. This truncation, which appears to generate a quasi-ideal Bose gas character to the composite 
boson system, mirrors the behavior of the well established ground state. This approach might not be suitable for other composite 
boson scenarios, which do not evolve directly from the BCS phase. Nevertheless, these crossover schemes provide a useful 
way of learning about composite boson systems in general. Moreover, they provide valuable insights about how to extend BCS 
theory slightly, without abandoning it altogether, and in this way to address a large class of short coherence length, but otherwise 
conventional, superconductors. 
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APPENDIX A: EVALUATION OF THE VERTEX CORRECTIONS 

In this appendix we demonstrate an explicit cancellation between the Maki- Thompson (MT) and Aslamazov-Larkin (AL) 
diagrams of Fig. 3. In this way we prove that the contribution to the vertex correction 5K from the superconducting order 
parameter is given by the Maki-Thompson diagram, and the pseudogap contribution 5Kpg comes from Maki-Thompson (MT) 
and Aslamazov-Larkin (AL) diagrams. It is easy to demonstrate a cancellation between the MT diagram and the AL diagrams, 
which will greatly simplify the calculations. In general, we have 

o 2 

5K^^g{K,K- Q)q, = -{MTU + Y.'^,,{P)Go{P k)^^^^^^ • q, (Al) 

where {MT)pg refers to the MT diagram contribution, and Tpg{Q ^ 0) is the T matrix or pair propagator. 
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To prove this cancellation, we notice that the vertex corrections in the (four-)current-current correlation functions can be 
obtained from proper vertex insertions in the single particle Green's functions in the self-energy diagram. In the pairing approx- 
imation {GqG scheme) we have 



^pg{K) ^Y.'^^aiK + L)G^{L)ip\^_ 



.L)/2' 



(A2) 



where L is the four-momentum of the fermion loop, this procedure leads to one Maki-Thompson diagram and two Aslamazov- 
Larkin diagrams. 



Obviously, L in the Eq. ( A2) is a dummy variable so that its variation does not change Tj{K), namely, 
= E [^P9{{K + L + AL)G,{L + AL)^\k-l-i^L)I2 " ^skiK + L)G^{V)^\k_^^i^ 

L 

^Y,^[r^^{K + L + AL) ~ Tpg{K + L)] Go{L + AL)^^^^^^^^^^ 



+%g{K + L) [Go(i + AL) - Go{L)] vIk-l-^l)/2 + %a{K + L)Go{L) 
Using G{K)G-^{K) ^ 1, we obtain 



f(K-L-AL)/2 f(K-L)/2 



} (A3) 



G{K + AK) - G{K) = -G{K)[G-^{K + AK) - G-^{K)]G{K + AK) 
= -G{K)K^{K + AK, K)G{K)AK^' , 

where G-^{K + AK) - G-^{K) w K^,{K + AK, K)AKt' is the full vertex. Similarly, we have 

G„{K + AK) - G„{K) = -Go{K)[Go\K + AK) - G^\K)]Go{K + AK) 
= -Go{K)X^{K + AK, K)GoiK)AK'^ , 



(A4a) 



(A4b) 



where G^^iK + AK) - Gq^{K) k \^{K + AK, K)AK'' is the bai-e vertex, and M'{K + AK, K) = (1, Vkek+Ak/2)- 
Equations ( A.4 ) correspond to the vertex insertions diagrammatically along the full and bare Green's functions, respectively. 
Using Tpg{K + L)= g/[l + gx{K + L)], we obtain 



Tpg{K + L + AL)- Tpg{K + L) = -Tpg{K + L + AL)[x{K + L + AL) - x{K + L)]Tpg{K + L 
Writing xiK + L) ^ El' GiL')Go{K + L- L')^i'-(K+L)/2' we have 

X{K + L + AL) -x{K + L) = Y, G{L') [[G„{K + L - L' + AL) - Go{K + L- L')] vI,_^k+l^ 



(A5) 



-AL)/2 



2 2 

^L'-{K+L+AL)/2 ~ ^L'-(K+L)/2 



]] 



+Go{K + L-L') 
On the other hand, writing x{K + L) = Y,l' G{K + L- L')Go(i')¥'^x+L)/2-L" we get 

X{K + L + AL) - x{K + L)= J2{[G{K + L - L' + AL) - G{K + L - L')] Go(L')¥'^k+l+al)/2-l' 



(A6) 



+GiL')GoiK + L - L') 



Vl'-{K+L-AL)/2 Vl'-(K+L)/2 



} 



(A7) 



Combining Eq. ( A.6) and Eq. (A7), we obtain to the first order of AL 



xiK + L + AL) -xiK + L) = ^Yl {^(^') ^^o{K + L - L' + AL) - G„{K + L- L')] ^1,_^k+l+al)/2 

L' 

+ [G{K + L-L' + AL)-G{K + L- L')] Go(i')¥'L-(K+L+AL)/2} , 



(A8) 



where we have assumed in general Lp]^ = ip'^f.. Substituting Eq. (A8) and Eq. (A5) into Eq. (A.3), we obtain 
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= - 



i ^rj,g(if + L + AL)Tpg{K + L) {g(L') [Go{K + L - L' + AL) - GoiK + L - L')] ^l._ 

LL' 



{K+L+AL)/2 



- [G{K + L-L' + AL)- G{K + L- L')] Go{L')ipl,_ 

- ^ Tpg{K + L) [Go(L + AL) - Go (L)] ^\k-l-al)Ii 



(K+L+\L)/2 



^Go{L + AL)^l 



K-L-AL)/2 



Y^Tp,{K + L)G^{L) 



2 2 

¥'(A'-L-AL)/2 ~ V(K-L)/2 



(A9) 



Comparing this with the analytical expressions corresponding to the diagrams in Fig. 3, it is easy to identify the first two terms 
as the two AL diagrams (which we denote by ALi and AL2) and the third one with the MT diagram for the pseudogap vertex 
corrections. Therefore, 



1 



[{ALi) + {AL2)] + {MT)pg + J2%giK + L)GoiL) 



2 2 

'^{K-L-AL)/2 ~ V(K-L)/2 



= 



(AlO) 



Finally, we have 



6Ki^pg{K,K- AL)AL^ = (ALi) + (AL2) + iMT)pg 

= -{MT)pg -2Y,TpaiK + L)Go{L) 



dip 



(K~L-AL)/2 

~dL 



■AL 



(All) 



Changing variables K + L ^ P, AL -^ Q leads to Eq. ( A.1 ) 



The two contributions which enter Eq. (pTh result from adding the superconducting gap and pseudogap terms which are given 
respectively by 



SAsc{K, K-Q) = -A2^^k(^k-qGo(-X)Go(g - K)\{Q - K, -K) 



(A 12a) 



and 



M^g(i^, K-Q) = -Y, Tp3(F)(^k-p/2<^k-q-p/2Go(F - K)Gq{P + Q- K)\^^{P + Q-K,P-K) 



+ Y.Tpg{P)Go{P - K) 



dvl 



k-p/2-q/2 

dku. 



(A 12b) 



APPENDIX B: FULL EXPRESSIONS FOR THE CORRELATION FUNCTIONS P, Po AND Poo 

It is useful here to write down the component contributions to the different correlation functions in the electromagnetic 
response. After adding the superconducting and pseudogap contributions one finds for the current-current correlation function 



n 
m 



P = 2E 



El 



1 - 2f{Ek) 



2E^ 



+ f'iEi.) 



^|(V6k)(Vek) - \{^el){W^l) 



-2E/'(i?k) 



n^ 



Y.^nE. 



f72 - (q • VSk)^ 



(Vek)(Vek) 



k ^i 



f^2 _ (q . VSk)2 



^^(Vek)(Vek)--(Ve2)(V^2) 



(Bl) 



and for the current-density correlation function 



Po = -2f7E--^/'(£;, 



qV£;k 



^w 



172 - (q . v£;k) 



rVek: 



(B2) 



and finally for the density-density correlation function 
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Pi 



00 



-2E 



AM 



1 - 2/(£;k) 



2-Bk 



/'(i?k) 



k i?k)2 r!2 - (q . Vi?k)2 



(B3) 



In deriving the first of these we have integrated by parts to evaluate 

Vek + VS(ii:) 

k 

'i-2f[E^) 1 r .,^ . ,^ . 1 



n 
m 



K k 

A2 



%^^ 



2E^ 



f'iE^) 



ifiiiVe^) ■ (Vek) - ^(Ve^) . (Vipl) 



2E/'(^k)(Vek)-(Vek). (B4) 



These expressions are then used to evaluate Eq. (p8[) in the text. 
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FIG. 1. Temperature dependence of the excitation gap (A'^)^ " and the order parameter (A) (normalized at T = 0) for (a) weak coupling 
BCS, (b) moderate coupling, and (c) strong coupling. The dotted lines represent the difference of these two energy scales, corresponding to 
the pseudogap parameter Apg. A strong pseudogap develops as the coupling strength increases. 
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FIG. 2. Dependence of the cutoff momentum A, (where the pair excitation spectrum crosses into the particle-particle continuum), on the 
coupling strength on a quasi-2D lattice. Here 4i is the half bandwidth. Because of strong damping of the pair excitations, they are irrelevant at 
low g and BCS theory is valid in the shaded region. Outside this regime, pair excitations become important. 
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FIG. 3. Diagramatic representation of (a) the polarization bubble, and (b) the vertex function used to compute the electrodynamic response 
functions. Here the wavy lines represent T and it should be noted that the thin and thick lines correspond to Go and G respectively. The total 
vertex correction is given by the sum of the Maki Thompson (MT) and two Aslamazov-Larkin (ALi and AL2) diagrams. 



20 



^ 

^ 







ycxx 



k„/k. 



10 



g/Qc 

FIG. 4. AB mode velocity cjvp as a function of the coupling strength (main figure) for various densities characterized by ko/kp in 3D 
jellium. Plotted in the inset is the large g asymptote Coc/vf, versus ko/kp, which varies as {kp/koY''^, as expected. 
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FIG. 5. (a) Normalized AB mode velocity, c/6ta, on a 3D lattice with an s-wave pairing interaction for various densities as a function of 
g, and (b) the large g limit for c/Qta as a function of density n for fixed —g/Qt — 20. (Here 6t is the half bandwidth). The dashed line in (b) 
shows a fit to the expected low density dependence n^' . Plotted in the inset is the velocity of the phase and density coupled collective mode 
v^/6ta with the particle-hole channel treated at the RPA level, for the same n as in (a). 
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FIG. 6. Temperature dependence of the real (Re c) and imaginary (Im c) parts of the AB mode velocity for moderate coupling (solid lines) 
and weak coupling BCS (dashed lines) in 3D jellium with fco = 4fcF . The mode is highly damped as Tc is approached. Because of inaccuracy 
due to neglected T- dependent amplitude effects, the solid curves are cut off slightly below Te- 
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